


JHU-TIPAC 97003 

DFTT 15/97 

KIAS-P97002 



Finite temperature effects on the neutrino decoupling 

in the early Universe 



N. Fornengo^'*" *, C.W. Kim'"''' \ and J. Song*" ^ 

"r-^' ■ ^ Department of Physics and Astronomy, The Johns Hopkins University, 

0\ '• Baltimore, Maryland 21218, USA. 

0\ 



Dipartimento di Fisica Teorica, Universitd di Torino 
and Istituto Nazionale di Fisica Nucleare, Sezione di Torino 
pLn ! via P. Giuria 1, 10125 Torino, Italy 

cn ! 

'^ School of Physics, Korean Institute for Advanced Study 
Seoul 130-012, Korea 
^ ■ (February 1997) 

(N 

en 

g : Abstract 

o\ 

■^ I Leading finite temperature effects on the neutrino decoupling temperature 

o ' in the early Universe have been studied. We have incorporated modifications 

jis ■ of the dispersion relation and the phase space distribution due to the pres- 

et . ence of particles in the heat bath at temperature around MeV. Since both 

'"^ ' the expansion rate of the Universe and the interaction rate of a neutrino are 

reduced by finite temperature effects, it is necessary to calculate thermal cor- 
p\ ' rections as precisely as possible in order to find the net effect on the neutrino 

j^ ■ decoupling temperature. We have performed such a calculation by using the 

finite temperature field theory. It has been shown that the finite temperature 
effects increase the neutrino decoupling temperature by 4.4%, the largest con- 
tribution coming from the modification of the phase space due to the thermal 
bath. 
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I. INTRODUCTION 

The standard hot big bang model appears to be a rehable description of the evolution of 
the early Universe, one of the most remarkable successes being the prediction of the present 
abundance of light chemical elements from the primordial nucleosynthesis calculations [|l|. 
The Universe is usually described as a hot, dilute gas of particles in nearly thermodynam- 
ical equilibrium [§,Q. During the early epochs, the particle species in the thermal bath 
underwent departure from the equilibrium one after another: one of the typical departures 
was that of neutrinos, when the temperature of the Universe was about 1 MeV 0. The 
neutrino decoupling has important (indirect) effects on the evolution of the Universe, since 
it happened at the time close to the neutron-to-proton ratio {n/p) freeze-out temperature 
T ^ 0.7 MeV and to the photon reheating by e^e~ annihilation (T < irte). The synthesis 
of light elements depends sensitively on the n/p freeze-out abundance which is determined 
by the interplay between the weak interaction rates and the expansion rate of the Universe. 
Both rates are influenced by the neutrino decoupling temperature [0,0. Neutrinos which 
were decoupled early do not share the entropy transfer with electrons, positrons and pho- 
tons in the medium. As a consequence, their temperature Ty becomes slightly lower than 
that of the other particles in equilibrium. On the contrary, if neutrinos are not totally de- 
coupled when the entropy transfer begins, they can share part of the e^ entropy and their 
temperature would be higher. A small change AT^ modifies the statistical distribution of 
neutrinos, and in turn affects both the weak interaction rates which maintain the equilib- 
rium between neutrons and protons, and the expansion rate of the Universe which is due 
to the change of the neutrino contribution to the total energy density. The overall effect 
is to shift the n/p freeze-out temperature and hence the n/p abundance when the nucle- 
osynthesis begins. The effect on the present Helium abundance Y was estimated to be of 
the order /\Y/Y ^ —Q.l{ATy/Ty) (the change is with respect to the calculation in which 
neutrinos are not reheated by electrons) [Q. Therefore, a precise knowledge of the neutrino 
decoupling temperature is desirable to gain a confidence in the estimates of the primordial 
element abundance. 

The standard calculation of the neutrino decoupling temperature is based on the as- 
sumption that particles in the thermal bath behave like free particles. The interactions are 
only responsible for the thermodynamical equilibrium, but do not contribute to the energy 
density of the Universe. However, particles in the medium feel effective potentials due to 
the interactions with other particles, which modifies their dispersion relations or introduces 
effective mass for the particle. In addition to this dynamical effect, the phase space available 
for the interaction is necessarily modified by the statistical distribution of particles in the 
medium. 

Our purpose is to examine whether or not the finite temperature effects can actually 
lower the neutrino decoupling temperature, leading to possible changes in the nucleosyn- 
thesis prediction of the present abundance of light elements. We explicitly and consistently 
include in the calculations the thermal effective mass of a photon and an electron, and the 
thermalized phase space distribution in the cross section in the Born approximation. Higher 
order interactions and radiative corrections to the neutrino interactions are not considered 
here, since we aim at the leading order corrections by the thermal bath. Another medium 
effect, which is not considered here, is the absorption or the emission of photons in the bath. 



This effect turns out to be important in the nucleosynthesis calculations for the reactions 
which involve three body initial or final state (such as neutron decay and its inverse decay), 
because the additional photon involved modifies sizably the phase space of the reactions 
^. In the case of two body reactions responsible for the thermal equilibrium of neutrinos, 
however, this higher order effect is small compared with the corrections calculated in the 
present paper, the overall magnitude of which will be shown to be of the order of 15% on 
the interaction rates, leading to a 4% shift in the neutrino decoupling temperature. 

The plan of the paper is as follows. In Section ||, we discuss how to include finite 
temperature effects in the calculation. In particular, we will briefly discuss the calculation of 
effective mass of a photon and an electron in the framework of Finite Temperature Quantum 
Field Theory (FTQFT), with the special attention to the MeV temperature range which is 
relevant for the neutrino decoupling. Limitations of the previous calculations in applying to 



the present problem will also be discussed. In Section |T|, we evaluate flnite temperature 
effects on the expansion rate of the Universe and the interaction rate of a neutrino, and 
show their effects on the neutrino decoupling temperature. 



II. THERMAL EFFECTS 

In the early Universe where the particles are propagating in a thermal bath, rather than 
in the vacuum, their dynamics and interactions are modifled to some extent. The behavior of 
particles in a thermal bath is systematically described in the framework of FTQFT 0. There 
are two equivalent formulations of FTQFT: the imaginary-time and real-time formalisms. 
In the present paper, we adopt the real-time formalism where the Feynman rules for all 
the vertices are identical with those in the vacuum, and the presence of the thermal bath is 
taken into account by the modiflcation of the tree-level propagators of fermions and bosons 
as 
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where a is the gauge-flxing parameter. In Eqs.(|l|) and (0), u^ is the four-velocity of the 
medium {u^ = (1,0) in the rest frame of the medium) and ripB are deflned as 



nF,B{x) = 0{x)fF,B{x) + d{-x)fF,B{-x) 



(3) 



where 9{x) is the step function and Jf^b are, respectively, the Fermi-Dirac (FD) and Bose- 
Einstein (BE) distribution functions 
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In Eq.(^ the (+) and (— ) signs refer to fermions and bosons, respectively, and t] is the 
chemical potential. In this paper r] = is assumed for all the species. 

Corrections analogous to that of Eq.(0) should also be applied to the propagators of 
massive gauge bosons, which are exchanged in the weak interactions of neutrinos. Due to 
the presence of statistical distribution functions, however, finite temperature corrections to 
their vacuum propagators are exponentially suppressed at T ~ 1 MeV <^ Mw, Mz- This 
reflects the fact that the bath is too cold to excite those very massive degrees of freedom. 

In Eq.(^), T is the temperature of the heath bath measured in the rest frame of the 
fluid. The presence of the bath does not violate the Lorentz invariance of the system since 
appropriate definitions of temperature, as well as of all the thermodynamical variables, can 
be obtained in any reference frame by suitable transformation laws 0J^. Since the early 
Universe as a thermal bath is conveniently described in the rest frame of the fluid itself, i.e. 
in the comoving reference frame, temperature T has a direct and simple meaning. We will, 
therefore, adopt the comoving frame throughout the paper. 

In the FTQFT formalism, the effect of the bath on the dynamical evolution of particles 
is taken into account by modifications to their dispersion relations, which can be recast in 
the definition of effective mass for the particle. (This effect can be evaluated by calculating 
the self-energy of the particle in the heat bath, and will be briefly reviewed in Sections 
|11 A| and UB ). The change of mass of the particle modifies its contribution to the energy 



density of the Universe and therefore its expansion rate. At the same time, it modifies 
the interaction rate due to the change of dispersion relations in the cross sections and the 
distribution functions f{E). 

In addition, the phase space distribution is also influenced by the background: the pres- 
ence of the same particles in the surrounding medium as those produced in the interaction 
processes reduces (enhances) the production probability for fermions (bosons), respectively, 
according to the statistics of the particles. The final state density factors are modified as 
follows: 

- — ,., ^ > - — — [1 — fF(E')] for fermions (5) 

{27r)^2E' {2'Kf2E'^ Jt\ )\ \ ) 

(2;0^^(2;0^[' + ^^("'^] for bosons. (6) 

In summary, the influence of the medium on the particle evolution in the Universe is 
due to the temperature-dependent shifts in the dynamical mass of the particles and the 
temperature-dependent modification of the interactions between particles. Both modifica- 
tions will be included in the calculation of the neutrino decoupling temperature in Section 
[m| . We now turn to the discussion of thermal mass of the photon, electron and neutrino at 
T ~ MeV. 



A. Effective Mass of the Photon 

In the thermal bath at T ~ MeV, a photon propagates through a medium made of 
electrons, positrons and neutrinos. Its propagation is therefore influenced by the interactions 
with the e"*" and e~ . The effect of these interactions on the dynamical evolution of the photon 



is taken into account by calculating the self-energy of the photon in the e^ background. The 
one-loop self energy diagram gives 

n'^'^(A;) = nr(A;) + nf(fc), (7) 

where k'^ is the photon four momentum, U.Q'^lk) is the vacuum polarization tensor at T = 
and 



(8) 
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describes finite temperature corrections. In Eq.@, p'^ is the four-momentum of the electron 
in the loop and tuq denotes the electron mass in the vacuum. The separation of the self 
energy into two parts, one referring to the T = case and the other coming from the presence 
of the medium, is attributed to the separation of the electron propagator in Eq.(0). The 
function Il'^'^{k) is divergent and, as usual, is subject to the electric charge renormalization. 
The function Ui^lk) induces a finite shift in the photon propagator, generating effective 
mass for the photon. The vacuum polarization tensor can be decomposed as [0,^ 

n'''' = 7TTik,w)P'"' + TrL{k,w)Q'"' , (9) 

where P'^^ and Q'^'" are orthogonal projection operators (for their explicit form, see [[7[) and 
ttt and ttl are scalar functions given by 
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k 
tttCK w) = -\iTL{k, w) + ^g^'Ti^, . (10) 



In Eqs.(^ and (|TOD, w = k'^u^ and k = Jw"^ — k^^k^. In the comoving frame, w{= k^) is 

the energy of the photon and k{= \k\) denotes the magnitude of its three-momentum. The 
decomposition in Eq.(^ allows one to write the propagator as (in the Feynman gauge) 

A^"" = =^ . (11) 

In the comoving frame, the P'^'^ term describes the transverse modes of the photon field, 
while the Q'^'^ term is a linear combination of longitudinal and time-like modes. Equation 
(|Tl|) shows that the propagator of physical (transverse) photon modes has a pole at k"^ = ttt- 
This is interpreted as a thermal-generated dynamical mass. The pole k"^ = n^ in the non- 
transverse part of the propagator describes the Debye screening length of the photon in the 
e^ plasma 0,0. 

Since we are interested in the leading order 0{a) correction to the mass of the photon, 
the vacuum dispersion relations fc^ = for the photon and p"^ = m,^ for the electron can be 
used in the right hand sides in Eqs.(|^) and ([T0|). Therefore, the effective mass of the photon 
is 
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where /xq = ttiq/T and the function /i(/io) is defined as 



x^ 



h{fio) = / dx , fpix) . (13) 

Note that the photon effective mass in Eq. (|12D has been obtained without any restriction to 
T, hence vahd for all temperature. Equation (|l^) gives the correct limit m^ = for T = 
and is in agreement with the result [m'f^{T)] = (27raT^)/3 obtained in Ref. [0 in the limit 
of high temperature T ^ w,mQ. Obviously this limit does not apply to the early Universe 
where photons are in thermal equilibrium with the mean energy w ~ T. Moreover, since 



T ~ ttiq at the neutrino decoupling temperature, the limit jJ,o/x = mo/T ^ 1 in Ref. |T^ 
cannot be applied. 

The thermal mass of the photon is almost linearly dependent on temperature. The 
photon effective mass is 0.115 MeV at T = 1 MeV and 0.241 MeV at T = 2 MeV. Even 
though the photon is still relativistic, its contribution to the energy density of the Universe 
is substantially reduced. 

B. Effective Mass of the Electron 

As in the case of the photon, the dynamics of electrons in a thermal bath is also modified 
by the electromagnetic interactions with background photons and electrons themselves. The 
interactions with neutrinos are suppressed at T ~ 1 MeV, for they involve the exchange of 
heavy bosons W and Z. Therefore, the effect of the bath on the propagation of the electron 
is expressed by calculating the electron self-energy in the presence of the ambient e"*", e~ 
and 7. The electron self-energy at one-loop level becomes 

S(;^) = So(i^) + St(i^) , (14) 

where So(l^) is the electron self-energy for T = and its thermal correction is 



6{k^)fB{k') 6[{{p + kf - ml)fB{p' + e) 



Mi>) = -2vre^ / TTTTlh^d^ + f + ^o)7m] (15) 
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where p^ is the four-momentum of the electron and k^ denotes the momentum in the 
electron-photon loop. The one-loop self-energy modifies the electron propagator S{p) as 

S{p)-' =i>-m,~ S(i^) = [^ - mo - Mi>)] - T.T{i>) . (16) 

The standard So(|^) leads to the definition of physical mass of electrons and the wave function 
renormalization in the vacuum. The temperature-dependent self-energy S^-d^) produces a 



finite shift in the dispersion relation as []TT| 
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where the function /i(/io) is defined in Eq.(|T3D and 
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where dimensionless quantities have been defined as u = |p|/T, ep = Ep/T, x = \p — k\/T 
and e^ = {Ep — Ek)/T. The resuh of Eq.(|l^) is vahd for all temperature. It gives the 
correct result mf^ = mo at T = and agrees with the result of ||TTHT3| in the limiting cases 
T -C mo and T ^ mo- The last two terms in Eq.(]T7|) are relatively small at T <^ mo, 
for they are exponentially suppressed by the Fermi-Dirac or the Bose-Einstein distribution 
function. Around T ~ m,o, however, the third term becomes important and has to be taken 
into account (for instance, h{mo/T = 1) ~ 0.543). The last term Ja{p) in Eq.(p!7D intro- 
duces momentum dependence in the electron effective mass, which is important especially 
in the calculation of the neutrino cross sections with the momentum-dependent phase space 



distributions. This term has been calculated in the limit e^ ^ /io [|T2l. This approximation 
is only valid at high temperature (T^ ^^o); for the most significant contribution of e^ to 
the integral is from its mean values, i.e. e^; ~ 1. In the temperature range of our interest 
(T ~ mo) , this approximation cannot be adopted and we have to resort to a numerical cal- 
culation. This function is always negative, and monotonically increasing from the limiting 
value 



p^o Jo 
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up to limp^oo Ja{p) = 0. In order to see how significant the modifications to mo due to 
the Ja{p) term is, we plot in Fig.l the deviation of the m.'f calculation including the Ja{p) 
term from that neglecting Ja{p), as a function of temperature. The dashed line refers to the 
calculation for the case with p = (p), where (p) is the mean value of the momentum in the 
bath. The solid line referring to the case with p = corresponds to the largest contribution 
of the Ja(p) term, the maximum of which is 3.3 x 10^'^ at T ~ 2 MeV. The other thermal 
corrections in Eq. (|l7|) are always larger than that due to the Ja{p) term at least by one order 
of magnitude, as can be seen by comparing Fig.l to Fig.2 (in Fig.2, the relative correction 
|1 — mg''^(T)/mo| is plotted as a function of T). This shows that the Ja{p) term is negligible 
around the temperature T ~ MeV. We will therefore neglect this term in our analysis. This 
turns out to be a great simplification in the calculations, especially for the interaction rate, 
because the momentum dependence in the electron mass is avoided. 

Figure 2 shows that the thermal corrections to the electron mass at T ~ MeV is sizeable: 
at T = 1 MeV the electron mass increases by 4.1%, and at T = 2 MeV the correction is as 
large as 16%. 



C. Effective mass of the neutrino 

Neutrinos also acquire effective mass in the presence of a medium. In tlie temperature 
range of our interest, tlie contributions come from weak interactions witli tlie ambient elec- 
trons and positrons. The effective mass squared of the neutrino has been shown to be of the 
order of Gf{N — N), where Gp is the Fermi constant and N{N) is the number density of 
electrons(positrons) |T^. Since we expect N c::^ N, the modification to the neutrino effective 
mass at T ~ 1 MeV is practically absent. Therefore, we neglect finite temperature effects 
on the dispersion relation of neutrinos in the following. 

III. NEUTRINO DECOUPLING 

Neutrinos in the early Universe, like all other particles, are kept in thermodynamical 
equilibrium through their interactions with the particles in the heat bath. As long as its 
interaction rate F is larger than the expansion rate of the Universe H, neutrinos remain 
in thermal and chemical equilibrium. When F becomes smaller than H, due to the re- 
duced temperature of the heat bath and therefore the increased distance between particles, 
neutrinos start to depart from the equilibrium and subsequently evolve independently from 
the other species. Even though this decoupling is not a sharp event, we can define that it 
happens when 

F(T,) = H{Td) . (20) 

The temperature when Eq.(pOD holds is defined as the neutrino decoupling temperature T^. 
In the following Subsections we will discuss finite temperature effects on the calculation 
of Trf. Both F and H decrease as finite temperature effects are incorporated. A detailed 
calculation is therefore needed to determine whether or not the decoupling temperature 
actually increases or decreases, and to estimate the size of the effect. 

A. Expansion rate 

In the standard big bang model, the dynamical expansion of the early Universe is gov- 
erned by the Friedman equation H] 



H 



SttG 



-P 



1/2 



(21) 



where G is the Newton constant. Under the assumption that all the particles are in ther- 
modynamical equilibrium, the total energy density p is 

P = Y.9-^j^,E{p)fm, (22) 

i 

where Qi is the number of internal degrees of freedom of the particle species i. E and p are 
its energy and momentum, respectively, which are related by the usual dispersion relation 
E'^ = p^ + mf , where m, is the mass of the particle. 
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Equation (p^) leads to the following expression of H: 



H = ^^9'J'iT)T^ , (23) 

where g*{T) is the number of degrees of freedom at temperature T, defined by 



..CT-EdV^fd. fr^; . (24) 



where /ij = rrii/T. In Eg. (P^ , Tj is the temperature of the species i. For the species i in 
equilibrium, Tj is equal to the temperature T of the Universe, but after the decoupling, its 
temperature does not need to be the same as T (this is actually the case for neutrinos). Note 
that the dimensionless g*{T) is proportional to the energy density of the Universe in units 
of T^, which is effectively dominated by highly relativistic particles; since particles near the 
transition from the relativistic to the non-relativistic regime can also contribute to g*{T), 
they should be also taken into account in a precise calculation. When a particle species 
becomes non-relativistic, its contribution to the energy density and therefore to g*{T) is 
exponentially suppressed. Around the neutrino decoupling temperature (T ~ IMeV), the 
particles in thermal equilibrium are photons, electrons, positrons and Ve{i'e)- We will assume 
that muon and tau neutrinos had already been decoupled, so that they do not interact with 
the fe- However, they still contribute to the energy density. 

Now, let us discuss the effects of the thermal background. Since electrons, positrons 
and photons have effective mass larger than that in the T = situation, g*{T) and hence 
H{T) are reduced. Figure 3 shows the relative change [1 — gl'^ {T) / g^,{T)] as a function 
of temperature. Although the photon mass is a few tenth of MeV and the electron mass 
increases by 5 ~ 10% at the temperatures around 1 ~ 2 MeV, the effect on the number 
of degrees of freedom is small: g*{T) decreases by a factor of 0.1%. This is because in 
this temperature regime g*{T) is dominated by three neutrino species. Due to the relation 
between H{T) and g*{T) in Eq. (]23|) , the ensuing effect on the expansion rate is a factor of 



two smaller: AH/H ^ 0.5{Ag^,/g^,) ^ 5 x 10^^. As will be shown in Section IIIB , finite 
temperature modifications to F are two order of magnitude larger than that to the expansion 
rate, leaving very little influence of H on the calculation of T^. 

B. Interaction rate and neutrino decoupling temperature 

Neutrinos are kept in thermal equilibrium by the interactions with electrons and positrons 
in the heat bath. Focusing on the electron neutrino z/g, the interactions such as 

z/g + e^ < — > Ue + e^ (25) 

Ve + e+ < — > z/e + e+ , (26) 

are responsible for kinetic equilibrium, and annihilation and creation processes like 

z/e + Ue ^ — ^ e^ + e^ , (27) 
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maintain neutrinos in chemical equilibrium. Since the interactions which involve more than 
two particles are suppressed by additional powers of the small coupling constants, we will 
not consider them in the following. For definiteness, we will discuss in detail the interaction 
rate for the process of Eg . (|25|) , which has the largest cross section and therefore is dominant 
in the determination of the thermal equilibrium for v^- We will turn to the other processes 
at the end of this Section. 

The standard calculation of the interaction rate of a neutrino relies on a number of 
simplifying assumptions: (1) the electron is considered to be massless, i.e. mg = 0; (2) 
the energy distribution of the initial-state particles is neglected: two initial particles are 
considered to have mean-valued energies 

and the interaction is supposed to occur in the center of momentum frame (cos^ = —1). 
Under these assumptions, the interaction rate becomes 

T = N,{T)vM(y{{E),{uo)), (30) 

where vm = 2. The cross section a for the massless electron is 

where s = A{E){uj) and the definition of the constant A will be given later (see Eq. (^Tj) ) . 
Since the electron is assumed to be massless, we have 

7 4 

IE) = {uj)= , , T-3.15T (32) 

^ ' ^ ' 180C(3) ^ ^ 



and 

N,{T) = ^^geT' ~ 0.182 T^ , (33) 

where C is the Riemann zeta function and C(3) = 1.20. Consequently, F is 



/ T 

^2 T^5 ^ Q /lo XX 1 n-22 / 
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F ~ 0.385 AG'pT'' ~ 3.48 x lO""' ( t^ttt j • (34) 

Comparing Eq.(|33D with the expansion rate calculated for the massless electron {g*{T) = 
10.75), 

the decoupling temperature is 

10 



Td - 1.09 MeV . (36) 

The above standard analysis provides very rough estimates of the interaction rate and the 
decouphng temperature. For comparison, the calculation of the expansion rate for rrif. = 0, 
but taking into account the thermal distribution of the initial energies, gives 

r~0.1284AG'|r5 , (37) 



which is a factor of 3 smaller than the estimate in Eq.(|3J). The resulting decoupling tem- 
perature is 

Trf ^ 1.56 MeV , (38) 

which is 50% higher than the value given in Eg . (|36D . 

A detailed and precise calculation of the interaction rate is rather involved, but necessary 
to see whether or not finite temperature effects lower the neutrino decoupling temperature. 
For the following interaction 

iye{uJ,k) + e~{E,p) < — y Ve{^\k') + e'{E\p) , (39) 

the interaction rate is defined by 

r(l/ee" -^ Z/eC") = ^— ^ j ^^ 9efe{E) Qufui^) [(XT^m] , (40) 

where the number density N^{T) of neutrinos at temperature T is 

A^.(r)=^./^/.(a^). (41) 

That is, r is the thermal average of product of the cross section a and the MoUer velocity 



vm 15 . The cross section of the process is, including the thermal phase space. 



(42) 

and the Moller velocity is defined (for massless neutrinos) by 

where s = (p + k)"'{p + k)a is the total energy in the center of momentum frame of the 
colliding particles and its expression in the comoving frame is {6 is the angle between p and 
k) 

s = ml + 2Euo - 2pk cos 9 . (44) 
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The scattering amplitude for the process VeG~ — ^ Ve^~ comes from two diagrams with 
Z in the t-channel and W in the w-channel. Since mean energies of interacting particles 
are of the order of the temperature T ~ MeV(-C Mw^z)i we can express the averaged square 
amplitude in the low energy limit as 

\M\' = lQGl[{v + afip^KWk'^) + {v- af{p"^K){p'^k'^) - {v^ - a^fml{k'^k'^)] , (45) 

where Gp is the Fermi constant, v = gy + I, a = qa + i, and gy (gA) are the vector 
(axial- vector) weak coupling constant. 

As a first step, we neglect the thermal phase space in order to investigate the effect of 
electron thermal mass alone. With the standard phase space, the cross section is 



a 



^2 



s J 



1 4' 

A + B^ + C^ 



(46) 



where 



A = 4(a^ + av + v"^) 
B = 2{2a^ -av-v^ 
C=(a-v? . 



(47) 
(48) 
(49) 



Finally, F is to be obtained by taking a thermal average of avM- A simple expression of F 
in terms of Bessel functions, involving a one-dimensional numerical integration, has been 
obtained in the case of Maxwell-Boltzmann (MB) distribution for the initial particles |]15 



We do not adopt this approximation here, because it induces larger uncertainties than the 
finite temperature corrections in the present analysis. For instance, particles with energy 
E ~ (E) would have 4.7% larger interaction rate in the MB distribution than in the FD 
distribution. Therefore, we have to resort to an improved method []TB|. 

Since the only angular dependence comes from the relative angle 6 between p and k, we 
have 



d^p d^k = A-np^dp 2-nk'^dk dcos 9 



(50) 



with the kinematical limits (0, oo) for both p and k and |cos^| < 1. With the isotropic 
distribution functions, the annihilation rate is 



1 QeQu Gl r^ p^dp 
N^{T) Stt^ 127r io E 



fe{E) 



°° k'^dk 



UJ 



fu{co)I{p,k) 



(51) 



where 



1 /•! 

I{p, k) = - I dcos0 

2 j-i 



'i_!!^y 



A + 5^ 
s 



4' 



Y. ^i^i 



(52) 



where the constants a^ and the functions Jj are obtained from the angular integration as 
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and 









ai = ml{C + 3A- 3B) 






a^ = ml{B-3A) 




a3 = A 




04 = ml{3B -A-3C) 




as = ml{3C - B) 




as = -mfC , 


1 


= 2 


2 


= 2{ml + 2Euj) 


3 


= 2iml + 2Eiu)'' + -pW 


4 


2puj 


'ml + 2Euj + 2pu;' 
m2 + 2Euj - 2pcj 






1 
2puj 


1 1 




5 


ml + 2Euo - 2puj ml + 2Euo + 2puj 






1 
4pu; 


1 1 


6 


{ml + 2Euj - 2pujY 


(m2 + 2Eu} + 2pcj)2 



The above holds for all interaction rate as long as the involved electron has momentum- 
independent mass. The integrations over p and k in Eq.(pTD are performed numerically (all 
the numerical calculations are performed with a precision of 10~^ but only three significant 
decimals are shown). 

The interaction rate for m,e = rno, i.e. without any thermal correction, is plotted in Fig. 4 
(dotted line) as a function of temperature. The decoupling condition Eq.(^) is satisfied at 
the temperature 



Tdirrie 



mo) 



1.57 MeV . 



(53) 



The inclusion of the electron thermal mass has the effect to slightly reduce the interaction 
rate (see the dashed line in Fig. 4): around T ~ 1.6 MeV the reduction is 0.3%. As a 
consequence, the decoupling temperature increases by 0.1%, i.e. 



TJm. 



m, 



eS\ 



1.58 MeV 



(54) 



Let us now take into account the thermal phase space in the evaluation of the cross 
section, which is 



dn 



dV d^A;' 



(27r)32E' {2tt)^2uj' 



;[! - feiE'm - U^')]{27^Y6^'\p + k-p'-k'). 



(55) 



The three-momentum p can be integrated out by using the three-dimensional delta function 
to give 
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dn= / ^ k'^dk'd<j)k'dcosek'[l-fe{E')][l-Muj'mE + u-E'~u') , (56) 

where the angle 6k' is defined as 

(p + fc) ■ k' 

cos Ok' = -. -. . 57 

\p + k\\k'\ 

After a trivial integration over 0^/, we use the remaining one-dimensional delta function to 
perform the integration over cos 6^' 



--s^i^ 



d{E' + u'] "' 



d cos Oy 



[I - UE')][l - Uuj')]duj\ (58) 



with the constraint p+k=p' + k'. This gives 

dn = ^P{uj')duj' , (59) 

where we have defined 

Pico') = [l-ME'{co'mi-Uu')], (60) 

with 

E'{uj') = a-u' (61) 

a = E + uj (62) 

^ = |p + A?| = (p2 + A;2 + 2pk cos ^)^/2 . (63) 

The quantities a and ^ are fixed with the given initial state variables. The kinematical 
limits for the final state energy uj' are 

"»" = 2(« - A) C"') 

""-''^ - 2(a + ^) ' ^^^^ 

which can be obtained from the constraint | cos^^'l < 1. 

The relevant phase-space distribution P{uj') is a function of the energy u' and of the 
energies and momenta of initial particles. In the case of non-thermal phase space, we have 
P{uj') = 1. The function P{uj') represents the reduction of the phase-space due to thermal 
effects. Figure 5 shows as a function of u' P{lj',E = {E),u = (u)) with cos^ = at 
three different temperature T = 0.1, 1, 5 MeV. We observe that the reduction due to the 
thermal distribution function is sizeable over all the energy range of u'. Especially around 
the temperature T ~ MeV, the amount of the reduction is always greater than 8% for all 
energies. 

Unfortunately, the functional form of P{uj') prevents us from analytically performing the 
integration over the final state energy u' and the initial state angle 6, contrary to the previous 
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situation without the thermal phase space. In order to minimize numerical integrations, we 
approximate the effect of the thermal phase-space as a mean value effect 

4Euj{crvM)TH = ^ f P{uj')du'\M^\ 
an A J 

~ ^{Pico')) fdcu'lM'l = AEu{avMUP{uj')) , (66) 

mA J 

where the quantity Euj{avM)o is the same as the one calculated previously (see Eqs.(^) and 
(1^)). In Eq.(|66|) the mean value of the phase space is defined as 

(P(cu')) = R{T; E,uj) = - j\ cos 9- ^—^ T"'" P{uo')du' . (67) 

The reduction factor R{T\ E, uj) is a function of the initial state energies and the temperature. 
In Fig. 6, R(T;E,u!) is plotted against the temperature for E = {E) and u = {u). As 
the temperature increases, R decreases because the Fermi blocking at higher temperature 
obstructs scattering processes. For instance, R{{E), (u)) = 0.86 at T = 1 MeV. Note that 
the effect even at T = 0.1 MeV is about 8% increase, where the finite temperature effects 
due to the thermal mass are totally negligible. 

With the definition of the reduction factor R{T; E,uj), the interaction rate including all 
the thermal effects is 

1 QeOu Gp f^ p'^dp f^ k'^dk 



N^{T) 87r4 127r Jo E 



r ^fe{E) r -^M^) Hp, k)R{T; E, u) . (68) 

Jo E Jo oj 



The interaction rate of Eq.(p8D is plotted in Fig. 4 as a solid line. The global reduction of F, 
compared to the situation without the thermal phase space, is 13% for T ~ MeV. This is 
by far the most important effect due to the presence of the thermal bath, being more than 
a factor of 40 larger than the inclusion of the electron thermal mass alone in the calculation 
of F. We recall that the thermal mass affects the expansion rate much less, leading to 
AH/ H ~ 5 X 10"'^ in the MeV range of the temperature. The decoupling temperature 
obtained from F in Eq.(^) is, finally, 

Td{me = m^^, thermal phase space) = 1.643 MeV , (69) 

which is a 4.4% increase from the decoupling temperature without any thermal effect. We, 
therefore, conclude that the thermal bath has the effect to increase the neutrino decoupling 
temperature, assuring that neutrinos are totally decoupled at the time of e^ entropy release. 
For completeness, we conclude this Section by reporting the neutrino decoupling tem- 
perature taking into account all the reactions listed in Eqs.(p5|-|27D. Including the thermal 
mass and the thermal phase space as discussed above, the neutrino decoupling temperature 
increases by 4.4% into 

Td - 1.41 MeV , (70) 

as compared to Td — 1.35 MeV obtained without finite temperature effects. 
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IV. CONCLUSIONS 

We have studied leading finite temperature effects on the neutrino decouphng temper- 
ature in the early Universe. The major motivation is to investigate if finite temperature 
effects could actually lower T^ and eventually affect the nucleosynthesis calculations. Two 
major features of the finite temperature effects have been incorporated in the calculation: 
(1) the interactions among the particles in the medium affect their dispersion relations which 
are recast in the form of effective mass ; (2) the presence of the medium modifies the phase 
space distribution of the particles in the processes that determine the equilibrium. 

The effect of the inclusion of thermal mass into the expansion rate of the Universe H has 
been shown very small: AH/H ~ 5 x 10~^. This effect turns out to be negligible compared 
to the modification of the interaction rate of a neutrino due to finite temperature. The 
latter has been discussed in detail for the process z/g + e ^^ z/g + e. The incorporation of 
the electron thermal mass alone leads to a 0.1% increase of T^. When the thermal phase 
space in the Born approximation is also considered, the decoupling temperature further 
increases. The total thermal effect is an increase of T^ by 4.4%, and the actual value of 
the decoupling temperature is T^^ (thermal effects) = 1.41 MeV. In conclusion, it is still valid 
even in the presence of a heat bath that neutrinos are totally decoupled at the time of e^ 
entropy transfer. 
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FIGURE CAPTIONS 



FIG. 1 The deviation of the effective mass of the electron ml^ calculated by including the 
Ja{p) term from the one without the Ja{p) term as a function of T. The solid line 
refers to the electron momentum p = 0, the dashed line to the mean value p = (p). 

FIG. 2 The deviation of the thermal electron mass nfj^ {T) from itlq as a function of T. 

FIG. 3 The deviation of the number of degrees of freedom with finite temperature corrections 
gf^iT) from 5'*(T) as a function of T. 

FIG. 4 The expansion rate H (dash-dotted line) and the interaction rate F for the reaction 
z/g + e~ < — > Ve + e~ as functions of T, in units of 10~^^ MeV. The dotted line refers 
to F when the thermal corrections are neglected. The dashed line represents the case 
with the electron thermal mass. The solid line includes both the electron thermal mass 
and the thermal phase space. 

FIG. 5 The phase space reduction function P{uj') plotted as a function of the final-state 
neutrino energy uj' and calculated for the mean values E = (E) and w = (w) of the 
initial-state energies with cos^ = 0. The three curves refer to different values of T: 
T = 0.1 MeV (dotted line), T = 1 MeV (sohd line) and T = 5 MeV (dashed fine). 

FIG. 6 The thermal phase space reduction factor R{T; E, u) as a function of T, calculated for 
E = {E) andcj = (lo). 
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